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LIST  OF  SYMBOLS 


velocity 

undisturbed  stream  velocity  at  the  airfoil  position 
coordinates  parallel  and  normal  to  the  stream 
reference  chord  length 
local  stream  shear 


pitching  moment  about  the  airfoil  midchord  station 

angle  of  attack 

density 

airfoil  thickness-chord  ratio 

maximum  height  of  the  airfoil  camber  line 

transformation  length 

parameter  defining  nonuniform  shear 

reference  velocity  at  the  stream  plane  of  symmetry 

stream  function 

half  height  of  the  two-dimensional  slipstream 
see  Eq.  (10) 

see  Eq.  (11) 

velocities  in  the  boundary  layer  parallel  to  and  normal  to  the 
surface 

pressure 

kinematic  viscosity 
boundary -layer  thickness 


displacement  thickness 
nnomentum  thickness 
see  Fig.  14 
A^/S 

velocity  at  the  outer  edge  of  the  boundary  layer 


INTRODUCTION 


The  Cornell  Aeronautical  Laboratory  is  conducting  a  program  of 
theoretical  and  experimental  research  on  low-speed  aerodynamics  as 
applied  to  STOL  and  VTOL  aircraft.  The  objective  of  this  program  is  to 
re-examine  certain  aspects  of  classical  aerodynamic  information,  in  the 
light  of  low-speed  flight  requirements,  with  the  aim  of  seeking  aerodynamic 
processes  which  might  be  exploited  to  enhance  low-speed  performance. 

One  aspect  of  propeller-driven  aircraft  which  has  recently  received 
increasing  attention  is  the  existence  of  strong  gradients  of  longitudinal 
velocity,  or  shear,  in  the  propeller  slipstream.  This  slipstream  shear 
interacts  with  a  wing  surface  and  can  alter  the  wing  characteristics.  In 
theoretical  treatments  of  a  wing  interacting  with  a  propeller  slipstream, 
the  first  important  simplification  is  the  replacement  of  the  slipstream  with 
an  ideal  uniform  jet,  free  of  all  velocity  gradients.  The  application  of  these 
theories  requires  that  one  equate  the  actual  slipstream  to  an  effective  uni¬ 
form  jet.  One  method  employed  is  to  assume  the  uniform  jet  has  a 
momentum  flux  equal  to  the  average  in  the  propeller  slipstream.  These  and 
similar  procedures  are  well  founded  on  momentum  considerations;  however, 
the  implicit  assumption  is  that  the  flow  nonuniformity,  the  shear,  does  not 
influence  the  wing  characteristics. 

In  the  interests  of  clarifying  this  aspect  of  the  wing -propeller  slip¬ 
stream  problem,  a  theoretical  and  experimental  investigation  was  con¬ 
ducted  to  quantitatively  determine  the  effects  of  stream  shear  on  airfoil 
characteristics.  In  the  theoretical  study,  the  problem  was  reduced  to  its 
most  elementary  form  in  order  to  isolate  the  inviscid  effects  of  shear.  This 
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investigation  was  then  followed  by  systematic  experimental  research  to 
check  and  to  extend  these  results  to  the  more  practical  situation  of  a  wing 
interacting  with  an  actual  propeller  slipstream. 

This  experimental  research  provided  the  desired  verification  and  ex¬ 
tension  of  the  theoretical  results.  In  addition,  the  results  showed  that  the 
airfoil  maximum  lift  varied  markedly  when  the  airfoil  was  positioned  at 
various  locations  in  the  shear  stream.  The  change  in  airfoil  maximum  lift 
with  stream  shear  was  explored  in  detail  in  a  subsequent  experimental 
study.  These  latter  experiments  showed  that,  in  one  instance,  by  reposi¬ 
tioning  the  airfoil  in  the  shear  stream  through  a  distance  comparable  to  the 
airfoil  thickness,  the  airfoil  maximum  lift  was  doubled.  Further  inviscid 
and  viscid  theoretical  studies  were  made  to  isolate  and  explain  this  destall- 
ing  phenomenon. 

The  purpose  of  this  paper  is  to  present  both  experimental  data  on  air¬ 
foil  maximum  lift  in  a  two-dimensional  shear  stream  and  theoretical 
studies  devoted  to  explaining  these  results.  In  this  connection,  the  avail¬ 
able  inviscid  theories  for  airfoils  in  uniform  and  nonuniform  shear  flows 
will  be  briefly  reviewed  to  illustrate  the  inviscid  effects  of  stream  shear  on 
airfoil  characteristics.  The  experimental  data  showing  the  destalling  phe¬ 
nomena  will  then  be  discussed,  and  the  results  of  the  theoretical  studies 
will  be  presented  to  illustrate  the  effects  of  shear  on  airfoil  stall. 
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THEORIES  FOR  AIRFOILS  IN  SHEAR  FLOW 


The  problem  of  the  influence  of  stream  nonuniformities  or  stream 
shear  on  airfoil  characteristics  is,  of  course,  not  new.  This  problem 
has  existed  virtually  since  the  introduction  of  the  airplane,  and  sporadic 
attention  has  been  given  to  the  problem  since  then.  Concerted  attention 
to  the  problem  has  not  been  warranted  until  now,  since  the  emphasis  in 
propeller-driven  airplanes  has  been  on  propellers  which  perturb  the  flow 
slightly,  and  the  slipstream  shear  was  accordingly  small.  The  current 
interest  in  diminishing  aircraft  take-off  and  landing  distances  has  had 
two  effects;  first,  it  has  made  take-off  and  landing  considerations  of 
paramount  importance,  and  second,  it  has  emphasized  propellers  pro¬ 
ducing  large  flow  perturbations.  This  latter  fact  implies  that  the  slip¬ 
stream  will  possess  large  shear,  and  consequently  requires  that  the 

influence  of  shear  on  airfoil  characteristics  be  re-examined. 

Uniform  Shear 

The  effect  of  uniform  stream  shear  was  first  treated  by  von  Sanden^ 

in  1912.  In  this  work,  consideration  was  restricted  to  a  wedge  profile  in  a 

stream  with  a  linear  velocity  gradient,  i.  e.  uniform  shear,  and  a  numerical 

2 

solution  was  obtained.  This  same  problem  was  taken  up  again  by  Tsien  in 

1943,  in  connection  with  the  influence  of  ground  level  gradients  in  wind 

velocity  on  airfoil  characteristics.  In  his  investigation,  Tsien  considered 

a  symmetrical  two-dimensional  Joukowsky  airfoil  in  an  inviscid  stream  of 

indefinite  extent  with  uniform  shear,  and  obtained  an  exact  solution  for  the 

airfoil  forces.  The  interest  in  deflected  slipstream  configurations  with 

3 

highly  flapped  wings  led  to  an  extension  of  Tsien 's  theory  by  Sowyrda  in 
1958  to  the  case  of  a  cambered  airfoil  in  a  stream  with  uniform  shear. 


) 


Again,  an  exact  inviscid  solution  was  obtained,  valid  for  an’^  thickness, 
camber,  and  angle  of  attack.  The  model  considered  both  by  Tsien  and  by 
Sowyrda  is  illustrated  in  Fig.  1.  One  important  feature  of  this  model  is 
that  the  entire  system,  both  airfoil  and  stream,  is  two-dimensional  and 
moreover,  that  the  stream  extends  indefinitely  far  above  and  below  the  air¬ 
foil.  This  is  in  contrast  to  a  wing-propeller  combination  in  which  both  the 
wing  and  stream  are  three-dimensional,  and  the  extent  of  the  shear  region 
is  finite.  The  second  important  feature  is  that  the  stream  shear  is  con¬ 
stant,  so  that  the  Joukowsky  transformation  can  be  applied  to  the  governing 
equations  of  fluid  motion.  This  single  fact  makes  it  possible  to  obtain  an 
exact  potential  solution  to  the  problem. 

The  exact  solutions  for  the  lift,  pitching  moment,  and  pressure  dis¬ 
tribution  are  given  in  the  original  papers  by  Tsien  and  by  Sowyrda,  and  will 
not  be  repeated  here.  It  is  shown  in  those  papers  that  the  drag  of  the  air¬ 
foil  is  identically  zero  since  the  airfoil  is  a  two-dimensional  closed  body  in 
inviscid  flow,  and  is  subject  to  no  spanwise  velocity  variations. 

Considerable  insight  into  the  effects  of  shear  flow  on  airfoil  charac-  I 

teristics  can  be  obtained  by  restricting  the  solutions  to  the  case  of  small  . 

airfoil  thickness,  camber,  and  angle  of  attack.  This  is  done  for  illustrative  ^  I 

purposes  to  show  the  shear  effects.  The  lift  coefficient  and  moment  coeffi¬ 
cient  about  the  airfoil  midchord  station  are,  to  the  first  order  in  thickness 
and  camber 
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(2) 


where  a  positive  pitching  moment  is  a  stalling  moment,  T  is  the  ratio  of 
the  airfoil  maximum  thickness  to  the  chord,  ^  is  the  maximum  height  of 

the  airfoil  camber  line,  and  k  is  the  local  stream  shear  (Fig,  1). 

*  Ug  el  y 

When  the  shear,  ,  is  taken  as  zero,  it  can  be  seen  that  Eqs.  (1) 

and  (2)  reduce  to  the  usual  thick  airfoil  results  in  which  the  slope  of  the 
lift  and  moment  curves  are  increased  because  of  thickness.  For  non-zero 
uniform  shear,  it  can  be  seen  that  the  shear  couples  with  both  the  airfoil 
thickness  and  camber  to  result  in:  (1)  an  increase  in  the  slope  of  the  lift  and 
moment  curves,  much  like  a  thickness  effect,  (2)  an  overall  increase  in 
lift  at  all  angles  of  attack,  analogous  to  a  change  in  camber,  and  (3)  an 
overall  increase  in  pitching  moment  analogous  to  a  distorted  camber  distri¬ 
bution.  These  changes  in  the  airfoil  characteristics  are  proportional  to  the 
stream  shear  and  can  be  viewed  as  effectively  associated  with  a  distortion 
of  the  airfoil  shape;  that  is,  there  is  an  equivalent  airfoil  which  would 
exhibit  these  same  characteristics  in  uniform  flow.  This  equivalent  airfoil 
would  be  thicker,  have  increased  camber,  and  a  different  camber  distribu¬ 
tion. 

It  will  be  noted  in  Eqs.  (1)  and  (2)  that  for  small  thickness  and  camber, 

the  thickness  and  camber  interactions  with  the  stream  shear  are  simply 

3 

superimposed.  The  exact  solution  shows  that  in  addition  there  is  a  coup¬ 
ling  of  thickness  and  camber  terms  which  does  not  appear  in  the  linearized 
form.  In  some  cases  this  coupling  term  is  as  large  as  either  of  the  other 
two  terms.  To  illustrate  this  and  the  magnitude  of  the  inviscid  shear  effect 


5 


AI-1190-A-7 


on  airfoils,  consider  a  17%  thick  section  with  11%  camber  in  a  shear  flow 
defined  by  k  =  5  .  This  airfoil  crudely  corresponds  to  one  with  about  50“ 
flap  deflection,  and  the  shear  corresponds  to  about  the  steepest  velocity 
gradient  found  in  a  propeller  slipstream.  It  is  shown  in  Ref.  3  that,  at  zero 
angle  of  attack,  with  no  shear,  the  airfoil  lift  coefficient  is  1.  6  .  With 

shear  defined  by  k  =  5  ,  the  lift  coefficient  at  zero  angle  of  attack  is 

ixi  4.  9  .  The  lift  increment  due  purely  to  shear,  AC^  3.  3  ,  is  about 
double  that  in  a  uniform  stream  due  purely  to  camber.  Moreover,  each  of 
the  shear  contributions  due  to  thickness,  to  camber,  and  to  coupling  of 
thickness  and  camber  are  about  the  same. 


Nonuniform  Shear 

The  theory  for  airfoils  in  uniform  shear  flow  provides  considerable 

insight  into  the  characteristics  of  airfoils  in  nonuniform  flows.  However, 

the  elementary  theoretical  model  of  the  flow  bears  little  resemblance  to 

that  occurring  in  an  actual  propeller  slipstream;  in  particular,  the  theory 

postulates  uniform  shear  whereas  the  shear  in  a  slipstream  is  markedly 

nonuniform.  Some  understanding  of  the  influence  of  nonuniform  shear  flows 

on  airfoil  characteristics  can  be  obtained  from  the  two-dimensional  inves- 

4  5 

tigations  made  by  Jones  ’  In  Ref.  4  Jones  considers  a  thin,  cambered, 
Joukowsky  airfoil  in  a  stream  with  a  parabolic  velocity  distribution  given  by 

where  ,  and  are  defined  in  Fig.  Z. 

The  parameter,  y-  ,  fixes  the  nonuniform  shear,  and  H  is  the  usual 
transformation  length  given  approximately  by 


H 


where  <c  is  the  reference  chord  length  and  T  is  the  thickness  chord  ratio. 
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In  contrast  with  the  uniform  shear  flow  problem  just  described,  this 
undisturbed  flow  is  characterized  by  a  stream  vorticity  which  is  not  con¬ 
stant,  but  varies  between  stream  lines.  From  Helmholtz's  law  it  is  known 
that  the  vorticity  must  be  conserved  along  a  streamline,  thereby  implying 
that  the  disturbance  stream  function  for  a  profile  in  this  flow  does  not  sat¬ 
isfy  Laplace's  equation.  This  is  the  fundamental  difficulty  in  treating  air¬ 
foils  in  nonuniform  shear  flow  since  it  precludes  an  exact  solution  in  closed 
form. 


Jones'  approach  to  the  problem  was  to  assume  the  airfoil  thickness 
and  camber  were  small,  and  that  the  stream  nonuniformities  were  small; 


1 


1  .  His  approximate  solutions  for  the 


in  particular,  ^  ^ 

airfoil  lift  and  pitching  moment  about  the  midchord  station  follow  after ‘neg¬ 
lecting  powers  and  products  of  these  quantities,  and  are  expressed  as  fol¬ 
lows. 


=  7/2  ^  Zfrcc^h  r-fCaa^tl  O) 


,4, 

In  the  limit  of  a  uniform  stream,  these  results  reduce  to  the  usual 
thin  airfoil  result  with  the  first-order  thickness  correction.  Additionally, 
these  approximate  solutions  demonstrate  another  effect  associated  with  the 
stream  nonuniformities;  that  is,  the  slopes  of  both  the  lift  curve  and  the 
pitching  moment  curve  are  increased  for  positive  values  of  the  shear  param¬ 
eter,  ^  .  To  illustrate  the  magnitude  of  this  effect,  consider  a  15%  thick 

airfoil  located  at  the  plane  of  symmetry  of  a  shear  stream  defined  by  ^  ^  1/4 
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(see  Fig.  2).  This  value  of  ^  is  the  upper  limit  on  the  range  of  validity 
for  the  theory,  and  physically  corresponds  to  a  situation  in  which  the  vel¬ 
ocity  parabolically  increases  by  about  1%  in  a  vertical  distance  equal  to  half 
the  airfoil  thickness.  This  shear  is  modest  in  comparison  with  that  occur¬ 
ring  in  propeller  slipstreams.  However,  Eqs.  (3)  and  (4)  show  that  this 
shear  results  in  a  20%  increase  in  the  lift  curve  slope  and  a  10%  increase  in 
the  moment  curve  slope.  A  shear  more  typical  of  a  propeller  slipstream 
was  experimentally  investigated  in  previous  research”'  The  shear  stream 
used  in  that  research  was  defined  by  ^<4?  2  .  Considering  a  15%  thick  air¬ 
foil  at  the  stream  plane  of  symmetry,  the  stream  velocity  increases  by 
about  7%  in  a  vertical  distance  equal  to  half  the  airfoil  thickness.  The  data 
reported  in  Refs.  6  and  7  show  that  the  shear  increased  the  lift  curve  slope 
by  250%  over  that  obtained  in  a  constant  velocity  stream.  The  conclusion  is 
that  in  the  two-dimensional  case,  nonuniform  shear  typical  of  that  existing 
in  propeller  slipstreams  can  cause  profound  changes  in  the  airfoil  lift  and 
moment  characteristics. 

In  comparing  Jones'  results  for  an  airfoil  in  nonuniform  shear  with 
Sowyrda's  uniform  shear  results,  it  will  be  noted  that  the  former  does  not 
predict  any  interaction  of  the  shear  with  the  airfoil  thickness  and  camber  to 
cause  a  camber-like  increase  in  lift.  Intuitively  one  would  expect  an  inter¬ 
action  similar  to  that  predicted  for  uniform  shear  since  at  locations  re¬ 
moved  from  the  stream  plane  of  symmetry,  the  stream  shear  becomes  fairly 
uniform.  This  speculation  is  qualitatively  substantiated  by  a  subsequent 

5 

two-dimensional  investigation  by  Jones  of  the  forces  on  an  elliptic  cylinder 
in  a  stream  with  a  hyperbolic  distribution  of  velocities.  These  results  show 
that,  to  the  first  order  in  airfoil  thickness,  the  lift  increment  due  to  shear 
is  given  by 
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A  =  ^/z 


r 


(5) 


where  ~  '\j  shear  in  the  undisturbed  stream  at  the  airfoil  posi- 

tion.  In  addition,  the  experimental  results  reported  in  Refs.  6  and  7  show 
an  overall  increase  in  lift  in  nonuniform  shear.  In  the  original  theoretical 
development  for  airfoils,  the  airfoil  thickness  and  camber,  and  the  stream 
nonuniformities  are  taken  to  be  small  and  the  products  of  these  terms  are 
neglected.  These  neglected  products  contain  the  shear  interaction  with  the 

airfoil  thickness  and  camber,  and  they  do  not  appear  in  the  final  results, 

© 

Eqs.  (3)  and  (4). 

The  shear  interaction  with  the  airfoil  thickness  and  camber  are  im¬ 
portant  for  two  reasons.  First,  the  interaction  results  in  an  overall  in¬ 
crease  in  lift  which  can  be  substantial.  More  important,  the  shear  inter¬ 
action  causes  an  additional  pressure  distribution  on  the  airfoil  which  would 
then  influence  boundary  layer  separation.  The  potential  importance  of  this 
interaction  with  respect  to  the  present  problem  led  to  a  partial  extension  of 
Jones'  theory  to  determine  the  effect  of  the  shear-thickness  interaction 
on  the  pressure  distribution  of  a  symmetrical  airfoil.  This  extension  is 
completely  reported  in  Appendix  I  and  will  be  briefly  summarized  here. 

Appendix  I  follows  Jones'  development,  but  retains  terms  involving 
products  of  the  airfoil  thickness  and  the  shear  parameter,  ^  ,  to  obtain 

a  partial  solution.  It  is  noted  that  the  primary  concern  is  with  the  tangen¬ 
tial  velocity  distribution  over  the  airfoil,  which  is  comprised  of  a  remote 
stream  component  and  a  disturbance  component.  The  partial  solution  con¬ 
sists  of  the  evaluation  of  the  additional  higher  order  terms  in  the  remote 
stream  component  of  the  velocity. 
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As  a  check  on  the  significance  of  this  partial  solution  including  the 
second-order  terms,  the  results  have  been  applied  to  determine  the  lift 
increment  due  to  the  thickness -shear  interaction.  The  resulting  lift  incre- 


Comparing  this  with  the  lift  increment  given  in  Eq.  (1)  and  Eq.  (5), 

AC,  r  ,  it  c.„  be  ==e„  «.at  partial  aolution  yield,  about  70%  oi 

the  increment  predicted  by  uniform  shear  theory  and  other  nonuniform  shear 
theories.  In  addition,  Eq.  (8)  is  found  to  agree  fairly  well  with  the  exper¬ 
imental  data.  This  provides  some  justification  for  subsequently  using  the 
partial  solution  for  the  additional  velocity  distribution  in  predicting  air¬ 


foil  pressure  distribution. 

To  summarize,  the  existing  theoretical  treatments  for  airfoils  in  non- 

uniform  flows  consider  the  two-dimensional  problem  in  an  unbounded  field. 

In  general,  the  treatments  fall  into  two  categories:  exact  solutions  for  air- 

2  3 

foils  in  flows  with  uniform  shear  ’  and  approximate  solutions  for  airfoils 

4  5  8 

in  flows  with  symmetric  nonuniform  shear  ’  '  The  literature  treating 

2  3 

uniform  shear  flows  *  shows  that  the  stream  shear  interacts  with  the  air¬ 
foil  thickness  and  camber  to  produce  an  overall  increase  in  lift,  much  like 
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an  increased  effective  camber,  an  overall  increase  in  pitching  moment 

analogous  to  a  distorted  camber  distribution,  and  a  small  increase  in  lift 

and  moment  curve  slopes  analogous  to  an  increase  in  airfoil  thickness.  The 

4 

theory  for  thin  airfoils  in  nonuniform  shear  provides  an  approximate  solu¬ 
tion  restricted  to  small  values  of  stream  shear  to  show  that  the  first- 
order  effect  of  the  nonuniform  shear  is  to  increase  markedly  the  slope  of 
the  lift  and  moment  curves.  Since  this  solution  neglects  products  of  the 

airfoil  thickness  and  camber  with  the  stream  shear,  the  shear-thickness 

,  5 

interaction  is  not  included.  However,  further  work  shows  this  interac¬ 
tion  should  be  present  in  nonuniform  shear  flows,  and  that  the  magnitude 
of  the  interaction  should  be  the  same  as  in  uniform  shear.  This  conclusion 

7 

is  substantiated  by  experiments  . 
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EXPERIMENTS  ON  AIRFOIL  STALL  IN  SHEAR  FLOW 


The  emphasis  in  the  experimental  research  has  been  to  systematically 
investigate  the  effects  of  stream  shear  on  airfoil  characteristics  by  experi¬ 
mentally  duplicating  a  flow  model  employed  in  the  theoretical  studies  and 
then,  progressively  departing  from  this  flow  model,  to  approach  a  more 
practical  flow  situation,  such  as  the  flow  in  a  propeller  slipstream.  This 
procedure  provides  a  valid  check  on  the  basic  theory,  but,  more  important, 
allows  one  to  isolate  a  number  of  complicating  features  of  the  physical 
problem.  The  inherent  difficulty  in  this  approach  is  that  closely  controlled 
shear  flows  must  be  produced  in  a  wind  tunnel. 

Apparatus 

The  technique  employed  to  produce  shear  flows  in  a  wind  tunnel  is  to 
insert  nonuniform  screens  or  grids  upstream  of  the  test  section  to  selec¬ 
tively  introduce  losses  in  the  stream  and  thereby  produce  a  stream  of  non- 
uniform  total  pressure.  Since  the  screens  are  nonuniform  and  since  the 
velocity  must  be  continuous  through  the  screen,  a  nonuniform  static  pres¬ 
sure  gradient  is  produced  at  the  screen  to  cause  flow  curvature.  This 
curvature  damps  out  with  distance  from  the  screen  because  of  the  wall 
constraint,  but  it  precludes  the  possibility  of  designing  the  screen  by  simply 
extrapolating  the  desired  velocity  distribution  back  to  the  screen  to  deter¬ 
mine  the  required  local  screen  losses. 

A  method  for  designing  nonuniform  screens  has  been  reported  by  Owen 

Q 

and  Zienkiewicz^  which  is  applicable  for  screens  which  perturb  the  flow 
slightly  to  produce  a  stream  with  a  linear  velocity  gradient.  However,  the 
present  research  is  concerned  with  large  nonuniform  and  discontinuous  shear, 
and  it  was  necessary  to  extend  that  work  to  include  these  cases.  This  ex¬ 
tension  and  the  methods  for  constructing  the  screens  are  given  in  Ref.  7. 

The  particular  screen  used  in  the  present  investigation  and  the  velocity 

calibration  are  shown  in  Fig.  3.  The  distribution  of  stream  shear  is  shown 

in  Fig.  4.  This  screen  design  was  selected  to  produce  a  simulated  two- 

dimensional  propeller  slipstream,  and  the  shear  within  the  jet  is  fairly 

< 

typical  of  that  behind  an  actual  propeller.  The  screen  was  constructed 
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of  metal  rods  of  1/8"  to  1/2"  diameter  spaced  in  a  metal  frame.  The  view 
shown  is  the  upstream  side  of  the  screen.  During  the  flow  calibration  it  was 
found  that  there  was  spottiness  in  the  flow,  evidently  stemming  from  the 
rapid  local  expansion  of  the  flow  to  an  equilibrium  condition  before  local  mix¬ 
ing  could  occur.  This  difficulty  was  eliminated  by  bonding  a  fine-meshed 
metal  honeycomb  to  the  downstream  side  of  the  screen  to  locally  constrain 
the  flow  and  allow  mixing  to  take  place. 

Calibration  data  are  shown  in  Fig.  3  for  the  illustrated  screen  and  are 
seen  to  agree  well  with  the  theoretical  curve.  The  largest  discrepancy  occurs 
at  the  edge  of  the  jet  and  reflects  the  effects  of  viscous  mixing  between  the 
two  streams.  The  flow  external  to  the  jet  is  seen  to  be  somewhat  nonuniform. 
This  is  ascribed  to  difficulties  in  maintaining  fabrication  tolerance.  The 
screen  solidity  required  to  produce  this  flow  is  about  85%,  and  small  devia¬ 
tions  from  the  design  solidity  produce  marked  changes  in  the  local  screen 
losses.  Since  the  present  interests  center  on  the  shear  within  the  jet,  the 
small  nonuniformities  in  the  external  flow  are  not  important. 

The  airfoil  model  used  in  the  experiments.  Fig.  5,  was  a  17%  thick 
symmetrical  Joukowsky  airfoil  which  completely  spanned  the  wind  tunnel. 

In  anticipation  of  future  experiments  in  an  axisymmetric  slipstream,  the 
model  was  segmented,  and  the  center  section  instrumented  with  a  three- 
component  strain  gauge  balance.  The  model  was  assembled  with  gaps  on 
either  side  of  the  instrumented  section  to  prevent  balance  interference.  The 
entire  model  was  then  wrapped  with  0.  005"  thick  sheet  rubber  to  prevent 
flow  through  the  gaps,  and  the  balance  was  calibrated  with  the  sheet  rubber 
in  place.  The  balance  was  designed  to  measure  a  maximum  lift  of  about  3 
pounds.  During  the  calibration  it  was  demonstrated  that  the  entire  system 
including  the  rubber  sheet  was  linear  to  better  than  0.  5%,  and  that  if  these 
nonlinearities  were  included  in  the  calibration,  the  balance  would  resolve  a 
lift  force  of  about  0.  01  pounds. 

The  experiments  were  made  in  the  subsonic  leg^^  of  the  CAL,  One-Foot 
High-Speed  Wind  Tunnel^  This  leg  of  the  wind  tunnel  has  a  test  section 
with  a  cross-section  of  17"  x  21",  and  is  operated  as  a  closed-throat,  non¬ 
return  type  tunnel.  It  operates  at  atmospheric  stagnation  pressure  over  a 
speed  range  of  0  -  180  fps. 
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This  experimental  work  was  devoted  to  measuring  the  airfoil  lift,  drag, 
and  pitching  moment  at  various  locations  in  the  shear  stream,  and  making 
observations  of  the  airfoil  boundary  layer  using  oil  film  techniques.  The 
tests  were  made  at  airspeeds  of  about  100  mph,  and  the  Reynolds  number 
based  on  airfoil  chord  length  varied  from  3  x  10^  at  the  center  of  the  stream 
to  6.  75  X  10^  near  the  edge  of  the  slipstream. 

Force  Data 

The  lift  .nd  pitching  moment  data  are  presented  in  Figs.  6,  7,  and  8.  The 
drag  data  are  not  presented  since  it  has  been  established^  that  the  artificial 
stream  turbulence  influences  the  viscous  drag.  The  effects  of  artificial  tur¬ 
bulence  on  lift  will  be  discussed  subsequently.  The  lift  data.  Figs.  6  and  7, 
have  been  reduced  to  coefficient  form  by  using  both  the  undisturbed  local 
stream  velocity  and  the  average  jet  velocity.  The  former  presentation  offers 
a  good  measure  of  airfoil  relative  efficiency,  while  the  latter  provides  a 
direct  comparison  of  lift.  The  pitching  moment  data  are  reduced  to  co¬ 
efficient  form  using  the  undisturbed  local  stream  velocity.  A  scaled  sketch 
of  the  airfoil  and  the  stream  velocity  calibration  is  included  in  Figs.  6,  7,  and 
8  to  show  the  position  of  the  airfoil  in  the  stream. 

The  angle-of-attack  range  covered  in  these  experiments  was  not  suffi¬ 
cient  to  obtain  the  usual  indication  of  airfoil  stall,  that  is,  a  definitive  maxi¬ 
mum  in  the  lift  curve.  However,  the  data  obtained  showed,  in  most  cases, 
that  stall  was  imminent,  and  the  figure  of  merit  used  here  is  the  lift  at 
which  the  stall  process  began,  termed  the  usable  maximum  lift.  This  was 
identified  using  both  the  pitching  moment  data  and  the  drag  data.  Stall  be¬ 
gins  when  there  are  large  regions  of  separated  flow  on  the  airfoil  which  grow 
with  increasing  lift.  The  growth  of  these  separated  flow  regions  is  manifested 
by  a  moving  center  of  pressure  on  the  airfoil,  and  by  increased  drag.  The 
criterion  used  with  the  moment  data  to  establish  the  beginning  of  the  stall 
process  was  the  point  at  which  ~  0  •  quantitative  criterion 

could  be  applied  to  the  drag  data,  and  a  sudden  increase  in  drag  was  inter¬ 
preted  as  the  beginning  of  stall.  The  lifts  at  which  these  two  measurements 
indicate  approaching  stall  are  shown  in  Figs.  6  and  7,  and  these  two  criteria 
are  seen  to  be  in  rather  good  agreement.  It  is  interesting  to  note  that  these 
criteria  also  agreed  well  with  the  observed  beginning  of  lift  unsteadiness,  a 
pilot  's  measure  of  approaching  stall. 
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The  most  important  property  exhibited  by  the  data  in  Figs.  6  and  7  is 
a  destalling  phenomenon  associated  with  airfoil  position  in  the  stream.  The 
data  show  that  with  the  airfoil  located  immediately  below  the  stream  plane 
of  symmetry  (-^  =  -  yg-  •  ~  g  ~  YE  ^  near  the  boundary  of  the  shear 

stream  ( -^  ~  ^  delayed  relative  to  the  observed  stall  at  other 

airfoil  positions.  The  magnitude  of  this  effect  is  evidenced  by  comparing 

y  1  1 

the  data  obtained  at  ~Y~  “  “  ft"  *  ^  ‘g  with  data  obtained  in  a  compar- 

able  uniform  slipstream'  '.  The  stream  shear  and  the  usable  maximum  lift 
coefficients  (based  on  average  stream  velocity)  are  summarized  in  Table  I. 


TABLE  1 

The  Maximum  Usable  Lift  Coefficients  in  Uniform  and  Nonuniform  Streams 


Stream 

Model  Position, 

Shear, 

Max.  Lift-Coefficient 

Uniform 

0 

0 

0.  70 

Nonuniform 

0 

0 

0.  75 

Nonuniform 

-1/8 

1 

o 

1.10 

Nonuniform 

+  1/8 

+  1.  78 

0.  48 

First,  in  comparing  the  usable  maximum  lift  measured  at  the  center  of 
the  uniform  and  nonuniform  streams,  it  is  seen  that  the  two  agree  closely. 
The  difference  in  lift  is  smaller  than  the  errors  in  the  corrections  applied 
to  account  for  the  differing  boundary  interference.  It  is  concluded  that,  for 
practical  purposes,  the  usable  maximum  lift  measured  in  the  center  of  the 
shear  stream  is  equivalent  to  that  measured  in  tlie  center  of  the  comparable 
uniform  jet. 

The  data  obtained  at  ~  +1/8  in  the  shear  stream  show  profound 

effects  of  stream  shear  on  usable  maximum  lift.  These  data  show  that  by 
’■'The  comparable  uniform  slipstream  was  free  of  velocity  gradients  with  a 
velocity  equal  to  the  average  velocity  in  the  nonuniform  stream  considered 
here.  The  resulting  changes  in  airfoil  lift  due  to  the  different  boundary 
effects  were  computed  and  applied  to  the  data  to  provide  a  valid  comparison. 
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moving  the  airfoil  below  the  stream  plane  of  symmetry  through  a  distance 
corresponding  to  half  the  airfoil  maximum  thickness,  the  usable  maximum 
lift  was  increased  by  about  50%.  Similarly,  in  moving  the  airfoil  above  the 
plane  of  symmetry,  the  lift  was  decreased  by  about  40%.  The  other  data  in 
Fig.  7  confirm  this  general  trend.  When  the  airfoil  is  in  the  vicinity  of 

aj  +3/8  to  1/2,  the  usable  maximum  lift  is  comparable  with  that  meas¬ 
ured  at  the  center  of  the  stream,  and  finally,  with  airfoil  near  the  slipstream 
boundary,  the  stall  is  evidently  delayed. 

It  is  interesting  to  compare  the  observed  maximum  lift  variations  with 
the  stream  shear  characteristics.  Fig.  4.  Below  the  center  plane,  where 
the  maximum  lift  is  augmented,  the  stream  shear  is  negative  and  the  shear 
derivative  is  positive.  At  the  center  of  the  stream,  where  the  maximum  lift 
is  unaltered,  the  shear  is  zero  but  its  derivative  is  positive.  Immediately 
above  the  centerplane,  where  the  maximum  lift  is  decreased,  both  the  shear 
and  its  derivative  are  positive.  At  +  3/8  to  1/2,  where  the  maximum 

lift  is  unaltered,  the  shear  is  a  maximum  but  its  derivative  is  near  zero. 

Near  the  slipstream  boundary  ( +  7/8),  where  stall  is  evidently  delayed, 
the  shear  is  positive  but  its  derivative  is  negative.  This  correlation  of 
maximum  lift  with  stream  characteristics  is  summarized  in  Table  II  and 
plotted  in  Fig.  9.  The  data  point  for  o  7/8  is  not  plotted  since  it  falls 
far  off  scale.  The  conclusion  reached  from  this  qualitative  type  of  correla¬ 
tion  is  that  the  observed  destalling  phenomenon  is  associated  with  both  the 

TABLE  II 


Correlation  of  Maximum  Lift  with  Stream  Shear 


Airfoil 

Position 

Shear, 

Shear 

Derivative,  kg 

Effects  on 
Maximum 

Lift 

Negative 

Positive 

Beneficial 

Negative 

0 

0 

Positive 

None 

0 

Positive 

Positive 

Detrimental 

Positive 

Positive 

~0 

None 

A,  0 

Positive 

Ne  gative 

Beneficial 

Negative 
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stream  shear  and  its  derivative,  and  it  appears  that  a  beneficial  destalling 
effect  is  obtained  when  the  product  of  the  stream  shear  and  its  derivative  is 
negative.  A  detrimental  effect  is  obtained  when  the  product  is  positive,  and 
sensibly  no  effect  is  observed  when  either  the  shear  or  its  derivative  is  zero. 

It  should  be  noted  that  apparently  this  same  destalling  phenomenon  was 

12 

encountered  by  Brenckman  in  his  wing-propeller  slipstream  experiment. 
His  experiment  consisted  of  measuring  the  spanwise  lift  distribution  on  a  , 
wing  of  infinite  span  intercepting  an  actual  propeller  slipstream.  The  wing 
was  located  so  that  it  was  nominally  centered  with  respect  to  the  slipstream, 
and  it  was  found  that  the  maximum  lift  coefficient  with  the  slipstream  was 
about  20%  greater  than  in  uniform  flow.  Brenckman  observed  this  destalling 
phenomenon  both  inside  and  outside  the  slipstream,  and  ascribed  the  phen¬ 
omenon  occurring  within  the  slipstream  to  be  associated  with  a  boundary 
layer  stabilization  mechanism.  His  experimental  model  precludes  the  pos¬ 
sibility  of  associating  this  mechanism  with  the  stream  shear  because  of  its 
three-dimensional  character. 


Oil  Film  Studies 

One  possible  cause  of  the  observed  destalling  phenomenon  in  the  present 
two-dimensional  experiments  is  the  occurrence  of  a  three-dimensional  sep¬ 
aration  process  with  accompanying  vortex  formations  which  could  energize 
the  flow  and  locally  delay  separation  on  the  instrumented  section.  As  a 
check  on  this  hypothesis,  the  complete  series  of  experiments  were  rerun 
using  oil  film  techniques  to  observe  the  boundary  layer  separation  process. 
The  aerodynamic  forces  were  monitored  during  these  experiments  to  in¬ 
sure  that  the  observed  phenomena  corresponded  to  the  previous  results. 
Typical  photographs  of  the  observed  separation  processes  are  presented  in 
Figs.  10,  11,  and  12. 

The  photographs  shown  in  Fig.  10  were  made  with  the  airfoil  at  the 
•  -3ll6  position.  The  view  is  looking  down  on  the  upper  surface  from 
a  point  about  one  chord  length  ahead  of  the  wing,  and  the  flow  is  from  the 
bottom  to  the  top  of  the  picture.  Early  in  the  original  research  ,  it  was 
found  that  with  the  wing  spanning  the  wind  tunnel,  the  wing  at  high  angles  of 
attack  caused  wall  boundary  layer  separation  which  resulted  in  unusual 


* 
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vortex  formations  on  the  wing.  To  prevent  this,  the  side  walls  were  per¬ 
forated,  and  suction  was  applied  to  bleed  off  the  wall  boundary  layer.  The 
wall  perforations  can  be  seen  in  the  photographs  in  Figs.  10,  11,  and  12. 

The  suction  and  distribution  of  the  perforations  were  experimentally  fixed  to 
provide  a  reasonably  two-dimensional  separation  pattern  . 

The  oil  film  technique  used  in  these  experiments  was  to  apply  oil 
colored  with  lampblack  at  the  airfoil  leading  edge  while  the  wind  tunnel  was 
running,  and  to  observe  the  patterns  after  an  equilibrium  pattern  had  been 
established.  A  separated  region  is  identified  as  a  region  where  the  oil  tends 
to  accumulate.  A  typical  sequence  of  this  operation  is  shown  in  Fig.  11. 
Photograph  6  shows  the  oil  pattern  shortly  after  application;  photograph  7 
shows  the  pattern  after  about  two  minutes;  and  photograph  8  shows  the  pat¬ 
tern  at  or  near  equilibrium.  In  monitoring  the  forces,  it  was  found  that  the 
oil  did  not  noticeably  influence  the  forces  at  any  time  during  the  procedure  , 
and  it  was  concluded  that  the  observed  formations  were  a  valid  indication  of 
boundary  layer  separation. 

The  important  conclusion  reached  from  Fig.  1  0  is  that  there  was  no 
substantial  separation  of  the  boundary  layer  up  to  an  angle  of  attack  of  12*. 
Several  isolated  separated  regions  can  be  noted  in  photograph  8  where  the 
oil  tends  to  accumulate.  At  oc  «  14.2*  ,  these  same  regions  are  apparent, 
and  at  oc  17.  3®,  the  oil  appears  to  accumulate  on  the  rear  third  of  the 
wing.  However,  there  are  no  regions  clearly  indicating  flow  separation. 

The  photographs  shown  in  Fig.  11  were  obtained  with  the  airfoil  posi¬ 
tioned  at  «  -1/8  ,  and  a  no-flow  photograph  is  included  for  comparison. 

At  =  8.  7®,  it  is  apparent  that  sufficient  time  was  not  allowed  to  reach  an 

equilibrium  pattern,  but  there  appear  to  be  no  regions  of  separated  flow. 

The  photographs  taken  at  oc  *  11.9°,  14.6°,  and  15.6°  show  five  or  six 
spots  of  oil  accumulation,  suggesting  local  flow  separations  that  do  not 
change  markedly.  At  o<  *  16.4°,  these  regions  appear  to  have  merged, 
♦Checks  were  made  at  small  angles  of  attack  to  determine  the  influence  of  * 
the  sidewall  suction  on  lift.  It  was  found  that  the  change  in  lift  coefficient 
due  to  suction  was  about  0.  01. 
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indicating  a  general  pattern  of  decreased  velocities  at  the  trailing  edge, 
but  still  no  definite  separation.  The  three-dimensional  patterns  appearing 
at  the  wing  extremities  are  caused  by  improperly  adjusted  bleed  air  through 
the  walls.  Subsequent  checks  showed  that  these  patterns  could  be  eliminated 
without  altering  the  general  oil  pattern.  The  photographs  in  Figs.  10  and  11, 
then,  qualitatively  substantiate  the  force  data  in  that  no  strong  flow  separa¬ 
tions  are  observed. 

The  photographs  shown  in  Fig.  12  were  obtained  with  the  airfoil 
positioned  at  *  +  1/8,  a  position  where  the  force  data  clearly  indicate 

the  beginning  of  stall  at  oi  jw  11®.  A  no-flow  photograph  is  shown  for  com¬ 
parison.  Again  at  cx  =  6.  3“  equilibrium  oil  patterns  had  not  been  fully 
established,  but  two  distinct  regions  of  separated  flow  are  apparent.  These 
are  probably  caused  by  some  three-dimensional  effect.  At  oc  =  9.  3®,  the 
separated  flow  is  evidenced  by  the  accumulation  of  oil  and  the  clear  regions 
near  the  trailing  edge  indicate  reversed  flow.  The  pattern  is  generally 
three-dimensional,  reflecting  the  sensitivity  of  separation  to  very  small 
spanwise  gradients,  but  a  fairly  clear  separation  line  is  evident.  At  oc  c  12.  3®,^ 
separation  has  moved  forward  to  the  midchord  position,  and  there  is  an  in¬ 
dication  of  leading  edge  separation  and  re-attachment.  At  oc  -  15.  3*  and 
17.  3®,  the  oil  patterns  show  a  leading  edge  separation  bubble  with  separa¬ 
tion  subsequently  occurring  at  the  quarter  chord.  A  pair  of  vortices  are 
evident,  and  cause  the  intricate  flow  patterns  in  the  separated  region.  It  is 
believed  that  the  position  of  these  vortices  is  influenced  by  wall  proximity, 

*  but  their  existence  does  not  appear  to  depend  on  the  walls.  This  is  sub¬ 
stantiated  by  the  photograph  at  oC  =  9.  3®.  Many  of  the  oil  accumulations 
indicated  weak  vortex  formations  normal  to  the  surface.  The  vortex  pair 
that  subsequently  formed  apparently  reflects  the  coalescing  and  strengthening 
of  these  distributed  vortices. 

To  summarize,  the  oil  film  studies  showed  that  no  unusual  three- 
dimensional  separation  patterns  were  present  on  the  wing  under  conditions 
corresponding  to  destalling,  and  that  a  normal  type  of  separation  was  present 
when  stall  was  observed.  The  conclusion  is  that  the  destalling  phenomenon  is 
associated  with  an  inviscid  or  viscid  interaction  with  the  nonuniform  stream. 

m 
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Thickness-Shear  Interaction 

The  final  item  to  be  noted  in  the  experimental  data  is  the  typical  shear 

effect  on  airfoil  lift.  The  data  in  Fig.  6  show  a  finite  lift  on  the  symmetrical 

airfoil  at  zero  angle  of  attack.  This  effect  was  discussed  in  the  previous 

section  in  connection  with  the  theoretical  treatment  of  airfoils  in  shear  flow. 

The  experimental  lift  increments  due  to  shear  are  tabulated  in  Table  III  and 

2  3 

are  compared  with  that  predicted  by  uniform  shear  theory  '  .  The  theo¬ 

retical  lift  increments  have  been  corrected  for  slipstream  boundary 

TABLE  III 


Comparison  of  Experimental  and  Theoretical  Lift  Increment  due  to  Shear 


Y/r 

AC^ 

p 

AC^(T5ien) 

ACL(correct  cd)^ 

-3/16 

-2.  60 

-0.  265 

-0.  695 

-0.  233 

-1/8 

-2.  07 

-0.  205 

-0.  554 

-0.  194 

-1/16 

-1,  15 

-0.  106 

-0.  308 

-0.  108 

0 

0 

0 

0 

0 

1/8 

1.  78 

0.  182 

0.  366 

0.  162 

3/8 

2.  77 

0.  236 

0.  569 

0.  198 

1/2 

2.  69 

0.  153 

0.  553 

0.  149 

7 

interference  and  these  corrected  values  are  also  tabulated.  It  can  be  seen 
that  the  corrected  theoretical  lift  increments  are  in  good  agreement  with  the 
experimental  data.  The  observed  discrepancies  are  believed  to  stem  from 
the  approximate  methods  used  in  accounting  for  the  boundary  proximity. 
Since  the  partial  extension  to  Jones'  theory'^  yields  substantially  the  same 
lift  increment  due  to  shear,  this  agreement  between  theory  and  experiment 
serves  to  verify  the  extension.  This  verification  is  important  since  the  ex¬ 
tension  of  Jones'  theory  will  be  used  subsequently  to  examine  the  airfoil 
pressure  distribution. 

There  are  a  number  of  potential  explanations  for  the  observed  destalling 
phenomenon  which  can  be  grouped  into  two  categories;  (1)  peculiarities 

^Corrected  for  slipstream  boundary  interference. 


associated  with  the  experiment,  and  (2)  effects  associated  with  the  stream 
shear  which  are  manifested  in  the  inviscid  and  viscous  flows  about  the  air¬ 
foil.  In  the  first  category,  there  are  the  effects  of  artificially  induced  tur¬ 
bulence  from  the  shear  screens  which  alter  the  airfoil  boundary  layer 
characteristics,  and  the  effects  of  the  slipstream  boundary  which  cause  a 
disturbance  flow  field  at  the  airfoil.  As  a  check  on  the  effects  of  artificially 
large  turbulence,  a  uniform  screen  made  of  1/16"  wire  with  a  solidity  cor¬ 
responding  to  the  average  in  the  shear  screen  was  mounted  in  the  wind 
tunnel  to  produce  comparable  turbulence  in  a  uniform  stream.  Experiments 
with  the  airfoil  model  showed  the  effect  of  this  turbulence  was  to  decrease 
the  drag  by  about  50%  and  to  increase  the  maximum  lift  by  about  8%.  The 
conclusion  is  that  the  turbulence  does  alter  the  maximum  lift  but  the  change 
is  small  in  comparison  with  the  observed  de stalling  effect. 

The  effect  of  slipstream  boundary  interference  was  previously  investi- 

7 

gated  by  testing  the  symmetrical  airfoil  at  various  vertical  positions  in  a 
uniform  two-dimensional  slipstream.  These  experiments  showed  that  the 
maximum  lift  varied  only  by  5%  when  the  airfoil  was  tested  at  positions  vary¬ 
ing  from  =  0  to  3/4.  The  conclusion  is  that  tiie  slipstream  boundary 

interference  is  not  important  in  determining  airfoil  maximum  lift. 

The  direct  effects  associated  with  stream  shear  appear  to  dominate  in 
influencing  the  maximum  lift,  that  is,  the  inviscid  effects  of  shear  that  alter 
the  airfoil  pressure  distribution,  and  the  viscous  effects  of  shear  interacting 
with  the  airfoil  boundary  layer.  These  two  topics  will  be  discussed  in  the 
next  section. 
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INFLUENCE  OF  STREAM  SHEAR  ON  BOUNDARY  LAYER  SEPARATION 


Within  the  framework  of  the  thin  boundary-layer  co.icept,  there  are  two 
effects  of  stream  shear  that  can  influence  boundary-layer  separation.  The 
first  is  inviscid  alterations  in  the  airfoil  pressure  distribution  stemming  from 
the  interaction  of  the  stream  shear  with  the  airfoil.  It  was  demonstrated 
previously  that  the  effect  of  uniform  shear  is  analogous  to  a  warped  airfoil 
in  a  uniform  stream,  and  the  additional  effect  of  nonuniform  shear  can  be 
viewed  as  an  effective  increase  in  airfoil  thickness.  Both  of  these  character¬ 
istics  imply  that  the  pressure  distribution  is  altered  in  shear  flow.  These 
effects  will  be  examined  in  this  Section  using  the  available  theories. 

The  second  effect  is  a  viscous  interaction  between  the  stream  shear 

and  the  boundary  layer  to  alter  the  boundary-layer  characteristics.  This 

13  1 6 

effect  has  been  considered  by  a  number  of  investigators  with  the  em¬ 

phasis  on  determining  the  changes  due  to  stream  shear  in  skin  friction,  con¬ 
vective  heat  transfer,  and  boundary-layer  thickness  with  no  external  pres¬ 
sure  gradient.  The  present  problem  is  concerned  with  conditions  influencing 

boundary-layer  separation,  and  hence  the  pressure  gradient  must  be  con- 

17 

sidered.  Boundary-layer  theory  for  bodies  in  a  uniform  stream  shows 
that  the  criterion  for  separation  depends  only  on  the  Reynolds  number  and 
the  pressure  gradient.  In  this  section  the  classical  laminar  boundary  layer 
development  will  be  generalized  to  the  case  of  a  stream  with  shear  to  deter¬ 
mine  the  effects  of  shear  on  the  criterion  for  separation. 

Pressure  Distribution 

The  experiments  described  were  made  in  a  stream  with  nonuniform 
shear,  and  consequently  consideration  will  be  restricted  to  the  two  theories 

4 

that  apply  for  nonuniform  shear  flow.  The  basic  theory  is  that  of  Jones 
which  is  restricted  to  thin  airfoils  and  to  streams  with  small  nonuniformities. 
The  restriction  to  small  stream  nonuniformities  precludes  the  calculation  of 
a  pressure  distribution  for  stream  conditions  corresponding  to  those  in  the 
experiment.  Similarly,  the  extension  of  Jones'  theory  (Appendix  I)  to 
retain  the  shear  interaction  with  the  airfoil  thickness  is  applicable  only  for 
small  stream  nonuniformities.  However,  a  theoretical  calculation  of  the 
pressure  distribution  will  demonstrate  the  trend  of  the  shear  effects  on 
pressure  distribution. 
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The  distribution  of  velocities  over  the  airfoil  given  by  Jones"*  has  been 
linearized  in  angle  of  attack,  neglecting  camber  effects,  as  follows: 

J  -  3 


whe  re 


U 


=  7-  jzl  ^  ^  ^  P  *■  '^9 


.  TliO 


X«/ 


-*  f(-^)  ~  >■  7(>&Z  « 

=  [2  +  ~  <3. 2  6  82  ^ 

[f 


oC 


«)  =  -T* 


/5,  =  0 


/t  =  [(8 


+  2  -  4 


-*-  8  'X^o  J 


'A 


(9) 


(10) 


'^o  ^  r/1.3 


"■-  '■'Tf  ^  ^ 


and  the  transformation  to  the  airfoil  coordinate  system  is 

'i. 

H 


'i-  ..  ^  .  O-tJp)'  (1^0 

-  *-i77r'^r 


2  Y-o 


(11) 
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The  additional  disturbance  velocity,  due  to  the  airfoil  thickness - 
stream  shear  interaction,  is  taken  from  Appendix  I  and  linearized  in  angle  of 
attack. 


MJ 


+  -2  2  f  -  2{o( 

-  oc  4^42  j^«>i  J  J  p 

-  ^2cc{-^J^  3fJ 


(12) 


The  formulas  given  in  Eqs.  (9)  and  (12)  have  been  used  to  calculate 

the  velocity  distribution  on  the  airfoil  used  in  the  experiments  (  f  *  0.  17). 

In  keeping  with  the  theoretical  restriction  to  small  nonuniformities,  the 

calculations  were  made  for  a  stream  defined  by  cr*  1/4  with  the  airfoil 

"f"  * 

positioned  at  yA  =  -  0.  2208.  These  locations  correspond  to  the  experi¬ 
mental  airfoil  positions,  f/r  =  t  I 

/8,  but  the  stream  shear  is  about  an  order 
of  magnitude  smaller  than  that  in  the  experiments.  The  calculations  were 
made  for  the  airfoil  at  a  lift  coefficient  of  0.  63  at  both  positions.  The 
velocity  distributions  are  shown  in  Fig.  13  and  are  compared  with  the  air- 
foil  velocity  distribution  at  the  same  lift  coefficient,  0.  63,  in  uniform  flow. 
First,  it  will  be  noted  that  the  effects  of  the  stream  shear  are,  in  general, 
to  cause  a  "saddle -backed"  distribution  in  that  the  peak  velocity  at  the  lead¬ 
ing  edge  is  increased  and  moved  forward,  and  a  velocity  peak  occurs  at 
the  trailing  edge.  In  addition,  by  changing  the  airfoil  position  from 
y/r  =  +  1/8  to  -1/8,  experimentally  shown  to  be  beneficial,  the  peak  vel¬ 
ocity  at  the  leading  edge  is  increased. 

This  behavior  contradicts  the  experimental  evidence  since  it  shows  the 
effect  of  shear  is  to  increase  the  adverse  velocity  gradient,  and  hence  would 
promote  separation.  The  cause  of  this  behavior  is  apparent  in  the  two  com¬ 
ponents  of  the  airfoil  velocity  distribution;  the  zero-thickness  component, 

Eq*  (9),  and  the  first-order  thickness  component,  Eq.  (12).  The  comparison 
in  Fig.  1  3  is  made  at  constant  lift  coefficient.  When  the  airfoil  is  located 
below  the  stream  plane  of  symmetry,  it  is  acted  upon  by  negative  shear. 


24 


AI-1190-A-7 


Referring  to  Eq.  (8),  the  negative  shear  causes  an  overall  decrease  in  lift, 
and  hence  (Eq.  12)  an  overall  decrease  in  the  airfoil  velocity  distribution. 

It  turns  out  that  this  decrease  is  nearly  uniform  over  the  airfoil  chord  with 
only  minor  changes  in  the  distribution. 

Since  the  comparison  is  at  a  constant  lift  coefficient,  the  decrease  in 
lift  due  to  the  negative  shear  requires  that  the  airfoil  be  at  a  larger  angle  of 
attack  when  located  below  the  stream  plane  of  symmetry.  This  increases 
the  magnitude  of  the  zero-thickness  velocity  component,  Eq.  (9).  It  can  be 
seen  that  with  the  airfoil  below  the  stream  plane  of  symmetry,  there  are  two 
effects  tending  to  cancel  each  other.  The  zero-thickness  effect  is  to  increase 
the  velocities  and  to  warp  the  distribution,  and  the  first-order  thickness 
effect  is  to  cause  an  overall  decrease  in  velocity.  Li  this  particular  instance, 
the  former  dominates  to  cause  increased  peaks  in  the  velocity  distribution. 

It  is  important  to  note  that  the  distributions  shown  in  Fig,  13  neces¬ 
sarily  correspond  to  a  stream  with  small  shear.  In  contrast,  the  experi¬ 
ments  were  made  in  a  stream  with  shear  that  was  nearly  an  order  of 
magnitude  larger.  It  appears  that  if  the  experimental  destalling  phenomenon 
is  to  be  explained  on  the  basis  of  airfoil  pressure  distribution,  it  will  be 
necessary  to  remove  the  theoretical  restriction  to  small  shear  and  to  retain 
higher  order  shear  effects.  This  conclusion  is  substantiated  by  the  experi¬ 
mental  data.  Referring  back  to  Table  II  and  Fig.  9,  it  was  noted  that  the 
qualitative  nature  of  the  destalling  phenomenon  suggested  that  the  dominant 
parameter  was  the  product  of  the  local  shear  and  the  local  derivative  of  the 
shear,  ^  J  ’  Ref'^rring  to  Eqs.  (9)  and  (12),  it  can  be  seen  that  there  arc 
no  terms  in  the  velocity  distribution  involving  this  parameter.  This  is  a 
consequence  of  the  approximations  used,  i.e.  ,  that  powers  and  products  of 
f .  f(^)  and  were  negligible.  Evidently  theoretical  treatments  in¬ 

cluding  the  effects  of  large  stream  shear  are  required  to  check  this  experi¬ 
mental  behavior. 

To  summarize,  the  approximate  solutions  for  airfoils  in  nonuniform 
shear  flows  show  two  effects  on  velocity  distribution  due  to  stream  shear. 

The  first,  a  zero-thickness  effect  indicated  by  the  average  uf  the  -y/r  s;  ^1/8 
curves  in  Fig.  13,  is  a  general  warping  of  the  distribution  to  cause  peaks  at 
the  leading  and  trailing  edges.  The  second  effect,  llu-  first-order  tliickness- 
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shear  interaction,  is  largely  fixed  by  the  local  stream  shear  both  in  sign  and 
magnitude.  Under  conditions  qualitatively  corresponding  to  those  where  the 
destalling  phenomenon  was  experimentally  observed,  the  sign  of  the  thickness 
effect  is  to  decrease  the  velocities  on  the  airfoil.  However,  this  effect  is 
more  than  counterbalanced  by  accompanying  changes  in  the  zero-thickness 
effect  so  that  the  net  change  in  the  velocity  distribution  is  detrimental.  The 
experimental  results  suggest  the  destalling  phenomenon  is  governed  by  the 
product  of  local  shear  and  the  derivative  of  local  shear,  a  parameter  which 
is  of  higher  order  than  those  considered  in  the  theory.  It  is  concluded  that 
the  destalling  effect  cannot  be  identified  with  the  existing  theories  which  are 
restricted  to  small-stream  shear,  and  that  solutions  valid  for  large-stream 
shear  are  required. 


The  Laminar  Boundary  Layer  with  External  Shear 

Following  the  thin  boundary-layer  approach,  it  is  possible  to  use 
theoretical  or  experimental  data  for  the  inviscid  pressure  distribution  to 
predict  the  separation  on  an  airfoil  if  an  appropriate  separation  criterion  is 

1  *7 

available.  Pohlhausen  '  has  established  a  criterion  for  bodies  in  uniform 
flow  using  momentum  integral  techniques,  but  the  applicability  of  this  cri¬ 
terion  to  the  present  problem  is  open  to  question  since  the  effects  of  stream 
shear  are  not  included.  The  momentum  integral  techniques  will  be  used 
here,  following  Pohlhausen,  to  determine  the  effects  of  nonuniform  stream 
shear. 


The  applicable  boundary-layer  equations  for  steady,  incompressible, 
two-dimensional  flow  are 


U. 


SU 


du  ^  J_ 

^  u  ^3  'V'  _ 

3  ^ 


(13) 

(14) 


where  U  and  -V”  are  the  velocity  components  along  and  normal  to  the 
surface,  and  y)C  and  are  coordinates  parallel  and  normal  to  the  sur¬ 

face,  and  i)  is  the  kinematic  viscosity. 

Integrating  Eq  .  (13)  from  the  wall  to  a  point,  jt  .  jvist  outside  the 
boundary  layer,  the  momentum  equation  is 
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(15) 


f  "1^ ^ ’’’I 

The  second  term  can  be  integrated  by  parts,  using  Eq.  (14),  to  arrive  at  the 
following  form 


pi 


ei'i^-UU) 


If  •'y 

Oa 


br  _ 


The  pressure  term  in  Eq.  (16)  can  be  evaluated  by  noting  that  the 
boundary-layer  equations  apply  for  the  outer  flow  and  that  the  pressure  just 
outside  the  boundary  layer  is 

Assuming  that  the  boundary  layer  is  thin  and  that  the  pressure  is  invariant 
through  the  boundary  layer,  the  pressure  term  in  Eq.  (16)  can  be  replaced 
by  Eq.  (17).  Integration  of  this  term  yields  the  following 

[ii^  -p  If  ‘‘r*  w)-om  Ij-  ‘'y  (18) 

o  • 

Substituting  Eq.  (18)  into  Eq.  (16),  the  resulting  equation  can  be  manipu¬ 
lated  further  to  arrive  at  the  following  form  of  the  momentum  integral 

It  can  be  seen  that  Eq.  (19)  is  similar  to  the  classical  momentum  integral, 
but  two  additional  terms  are  introduced  by  the  stream  nonuniformity;  the 
third  term  on  the  left  and  the  second  term  on  the  right.  The  magnitude  of 
each  of  these  terms  is  fixed  by  the  shear  in  the  external  stream.  If  the 
stream  shear  is  sniiill  in  comparison  with  the  shear  in  the  boundary  layer, 
these  two  additional  terms  can  be  neglected  in  comparison  with  the  remain¬ 
ing  terms.  In  the  present  experiments,  the  stream  shear  was  about  3%  of 
the  boundary- laye r  shear.  Hence,  for  present  purposes,  the  applicable 
momentum  integral  is  the  classical  relation. 


AI-  1  190-A-7 


The  first  two  integrals  in  Eq.  (20)  are  identical  to  the  classical  defin¬ 
itions  of  the  momentum  and  displacement  thickness,  thereby  suggesting  that 
these  definitions  would  be  applicable  in  nonuniform  flows.  The  displacement 
thickness  used  here  is  defined  in  the  usual  manner. 


n  «o 

I  (U-  u)  ot 


Similarly,  the  momentum  thickness  is  defined  as 


r®  *  r  “ 

j  U  =  u(U-u)<^y^ 

“*0  J/, 


The  present  problem  is  to  investigate  the  effects  of  a  stream  with  a 
parabolic  distribution  of  velocities.  In  this  connection,  consider  a  flow  with 
the  velocity  distribution  given  by 


This  is  the  same  external  flow  considered  by  Jones’,  but  the  coordinate 
is  now  referred  to  the  surface  of  the  body. 

The  flow  model  is  illustrated  in  Fig.  14,  and  since  the  external  inviscid 
flow  has  been  effectively  displaced  by  the  displacement  thickness,  the  equi¬ 
valent  inviscid  flow  over  the  surface  is 


Using  this  inviscid  velocity  distribution  in  Eqs.  (21)  and  (22),  the  displace¬ 
ment  and  momentum  thickness  are  given  by  the  following  relations 
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In  the  limit  of  uniform  external  flow  (  -  0),  these  relations  are  reduced  to 

the  usual  expressions  for  the  displacement  and  momentum  thickness.  By 
retaining  the  classical  definitions  of  displacement  and  momentum  thicknesses 
for  flows  with  external  stream  shear,  the  relations  have  been  considerably 
complicated  in  that  the  equations  are  cubic  and  quintic  in  the  thickness. 

Following  the  Pohlhausen  approach,  it  is  assumed  that  the  boundary- 
layer  profile  is  adequately  represented  by  a  quartic  series  expansion 


2  i 

(ar|  +  bq  +  ^rj  +  dq 


where  is  the  boundary-layer  thickness,  and  a,  b,  <c,  and  d,  are 

free  constants  to  be  determined  from  the  boundary  conditions.  The  appro¬ 
priate  boundary  conditions  are 
At'^sO:  ix  -  nf'  -  0 


At  =  S  : 


,  3*u  *  ..  30 


U  -  U  -  U. 


S  ^  6 -  S^ 

H 


^7  ^ 


9  u  ^30  _ 
3^"  h 


Applying  these  boundary  conditions  to  Eq.  (Z7)  the  velocity  profile  is 
determined. 


where 
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A  = 


$*_  -au 


All 


2  H  \  H  / 


Zr^-Zri^ri* 


~  l('l~  ^ 


H(q)  =  j(-5>ji-9fl^-^,j*) 


i(i)‘K-t^^i-'lV 


The  parameter,  /\,  ,  in  Eq.  (28)  is  the  familiar  Pohlhausen  shape 

factor  which  expresses  the  influence  of  pressure  gradient  on  boundary -layer 
development.  In  the  limit  of  a  uniform  stream,  ^  — >0  .  Eq.  (28)  reduces 
to  Pohlhausen's  result.  The  effects  of  the  stream  shear  are  introduced 
through  two  additional  shape  factors,  X.  and  ^  .  Noting  that  the  denom¬ 
inator  in  these  two  parameters  is  simply  / U,  ,  where  6/j  is  the 

velocity  at  the  outer  edge  of  the  boundary  layer,  it  can  be  seen  that  the  sum 
of  these  two  shape  factors  can  be  regarded  as  a  "shear  number"  or  "vorticity 

number".  That  is,  a  measure  of  the  shear  in  the  boundary  layer  is  Ug/S  • 

t-U  /  S  6  —  «r*  \  ^ 

The  shear  at  the  edge  of  the  boundary  layer  is  -£ — ^ _ j  .  Since 


A.  +  I  = 


\H  \  M  J 


it  follows  that  this  sum  is  the  ratio  of  the  free  stream  shear  to  the  boundary- 
layer  shear.  The  fact  that  these  two  parameters  do  not  appear  as  a  simple 
pair  implies  that  the  pertinent  stream  shear  is  not  that  at  the  edge  of  the 
boundary  layer,  but  rather  some  average  value. 

Using  the  results  of  Eqs.  (25),  (26),  and  (27),  the  momentum  and  dis¬ 
placement  thickness  can  be  written  down  immediately 
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-|- Ij- A- 25 •  s(-|^)  r£  (29) 

-f-|(l-X-2C)%2tO-A-2t)(-|-)*-|-  xO-A-21 

(30) 

_  37  A  _  A*  )2 

~  3  15  ~  9^5  907Z  *  105  ~  315  ^  *  Z835 

64  /3XA  /3  ,  533  ,  > 

315  ^  22GBO  3780  1890  ^ 

These  formulas  also  are  reduced  to  Pohlhausen's  results  in  the  limit  of  van¬ 
ishing  stream  shear  (  X  =  ^  =  O  ). 

The  important  consideration  in  the  present  problem  is  the  criterion  for 
separation.  Separation  occurs  when  the  shear  at  the  surface  vanishes.  Ap¬ 
plying  this  condition  to  Eq.  (28),  it  follows  that  separation  occurs  when  the 
following  criterion  is  met. 

A  ^  -12  ^  10  X  12  %  (32) 

It  is  possible  to  combine  the  relations  for  the  momentum  and  displace¬ 
ment  thicknesses,  Eqs.  (25)  and  (26),  with  the  momentum  integral,  Eq.  (20), 
and  cast  the  resulting  differential  equation  into  a  form  that  can  be  solved, 
using  Eq.  (31),  to  predict  separation.  This  procedure,  in  the  case  of  a 
uniform  external  flow,  is  feasible  because  of  the  elementary  form  of  the  re¬ 
lations  for  the  momentum  and  displacement  thicknesses  and  because  the 
velocity  profile  can  be  represented  by  a  single  parameter.  The  complicating 
feature  when  the  external  flow  is  nonuniform  is  that  the  expressions  for  tlie 
momentum  and  displacement  thicknesses  involve  powers  of  the  thicknesses, 

and  so  result  in  a  complicated  nonlinear  differential  equation.  Numerical 

# 

methods  could  be  used  to  solve  the  differential  equation;  however,  considerable 
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understanding  of  shear  phenomena  can  be  obtained  from  elementary 
considerations. 


The  approach  taken  here  is  to  use  a  linear  approximation  for  the  mo¬ 
mentum  and  displacement  thicknesses  in  an  approximate  calculation  follow- 
1  ft 

ing  that  of  Prandtl  .  This  approach  is  used,  not  to  arrive  at  a  quantitative 
criterion,  but  to  determine  the  qualitative  effects  of  stream  shear  on 
separation. 


Taking  the  displacement  and  momentum  thicknesses  to  be  small  in  com- 
parison  with  the  wing  chord,  Eqs.  (25)  and  (26),  to  the  first  order  in  ■ 


reduce  to 


''«o 

(U-u) 

•m) 


06  u  {U-  u) 


Inserting  these  relations  in  Eq.  (20),  the  momentum  equation  becomes 


(P^) 

/y--o 


Substituting  for  the  surface  shear  terms  using  Eqs.  (24)  and  (27)  and  re¬ 
arranging  terms,  Eq.  (35)  can  be  reduced  to  the  following  form 


U.Q  d6 


Following  Prandtl's  argument,  it  is  assumed  that  the  boundary  layer  is 
acted  upon  by  the  inviscid  pressure  distribution  up  to  a  point  which  is  close 
to  separation.  Starting  at  this  point,  it  is  assumed  that  the  pressure  gra¬ 
dient  is  such  that  the  shape  of  the  velocity  profile  remains  unchanged 
further  downstream;  that  is,  the  shape  factor,  ,  is  a  constant  from  this 

point  on.  By  definition. 


_  A  —L 

j  oLU 


and  after  differentiation 


e  de  A/0\*  ^6 
j  dK  "  z^s)  ruJ)* 
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where  the  primed  symbols  denoted  differentiation  with  respect  to  the 
x-coordinate.  Substituting  Eq.  (37)  into  Eq.  (36)  and  rearranging  terms, 
the  following  relation  is  obtained 


Using  Eq.  (38)  and  the  preceding  assumptions  as  to  the  boundary-layer 
characteristics  near  separation,  a  velocity  criterion  can  be  established 
for  preventing  separation.  Prandtl  does  this  by  assigning  a  value  to  the 
shape  factor,  /i-  =  -10,  which  is  near  the  value  for  separation.  The  equi¬ 
valent  procedure  here  is  to  take 


A  =  -  /O  ^  /o  A  +■  12  1 


(39) 


Substituting  this  value  into  Eq.  (38),  the  criterion  for  just  preventing 
separation  is 


2 

O/S 


± 

3 


/O-  I0\  -  ill 


(40) 


This  criterion  is  reduced  to  that  of  Prandtl  in  the  limit  of  vanishing  stream 
shear,  A  =  ^  =0.  The  significance  of  this  criterion  is  unchanged.  That 

is,  separation  can  be  delayed  if  the  second  derivative  of  the  velocity  is  suf¬ 
ficiently  large  and  is  positive.  Indeed,  a  prerequisite  for  preventing  sepa- 

^  II 

ration  is  that  Ug  be  positive.  K  it  is  positive  and  if  the  product  UgUg 
is  large  in  comparison  with  the  velocity  gradient  term,  )  >  separa¬ 

tion  can  be  delayed. 

Some  insight  into  the  significance  of  the  shear  in  altering  the  criterion 
for  separation  can  be  obtained  by  assigning  numerical  values  to  the  shape 
factors,  A.  and  ^  .  The  values  used  here  correspond  to  the  properties 

of  the  shear  stream  generated  in  the  experiments  with  the  airfoil  located  at 
y'/f'  =  ^  1/8.  For  purposes  of  comparison,  the  case  of  a  uniform  stream 

is  included.  The  displacement  and  momentum  thicknesses  were  calculated 
using  only  the  first-order  terms  in  Eqs.  (29)  and  (30)  and  are  summarized 
in  Table  IV. 
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TABLE  IV 


The  Effects  of  Stream  Shear  on  the  Separation  Criterion 


The  data  in  Table  IV  demonstrate  the  effect  of  shear  on  two  criteria  for 


separation.  One  criterion  is  that  previously  derived,  the  pressure  gra¬ 
dient,  y\.  ,  at  separation.  Eq.  (32)  and  Table  IV  show  that  the  parameter, 

+  A  ,  is  always  detrimental  and  promotes  separation.  However,  the  par¬ 
ameter,  ^  ,  can  take  on  a  negative  sign,  corresponding  to  the  experiments 

below  the  stream  plane  of  symmetry,  ^and  in  this  case  separation  is  delayed. 

In  contrast  the  criterion  on  77f~  shows  that  positive  A  and  F 

i^S) 

are  beneficial  in  that  a  wider  range  of  velocity  gradients  is  permissible. 

For  positive  A  and  negative  ^  ,  the  conditions  under  which  the  de¬ 

stalling  phenomenon  was  observed,  the  criterion  for  separation  is  nearly 
identical  to  that  for  a  uniform  stream. 


The  conclusion  is  that  stream  shear  can  alter  the  criteria  for  boundary- 
layer  separation,  but  under  conditions  corresponding  to  those  in  the  ex¬ 
periments,  these  alterations  are  negligibly  small  and  the  uniform  stream 
criteria  apply. 
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CONCLUDING  REMARKS 


In  summary,  the  two-dimensional  experiments  made  with  an  airfoil 
in  a  symmetric  nonuniform  shear  flow  disclosed  that  under  conditions  when 
the  product  of  the  local  shear  and  the  local  derivative  of  shear  was  nega¬ 
tive,  airfoil  stall  was  markedly  delayed.  When  this  product  was  positive, 
stall  was  prom.oted.  Subsequent  experiments  served  to  eliminate  extrane¬ 
ous  experimental  factors  as  the  cause  of  this  phenomenon  and,  by  elimina¬ 
tion,  to  isolate  the  cause  as  stemming  from  inviscid  alterations  in  the 
airfoil  pressure  distribution  and/or  changes  in  the  criterion  for  boundary- 
layer  separation. 

A  Pohlhausen  calculation  of  the  laminar  boundary  layer  with  a  non- 
uniform  external  stream  disclosed  that  the  external  shear  tends  to  alter  a 
criterion  for  separation.  However,  a  numerical  example  showed  that  in 
the  experiments,  this  effect  would  be  negligibly  small,  leading  to  the  con¬ 
clusion  that  shear  effects  to  alter  the  boundary-layer  characteristics  were 
not  the  source  of  the  destalling  phenomenon. 

The  available  approximate  airfoil  theories  were  applied  to  the  present 
problem  to  determine  the  effects  of  shear  on  the  distribution  of  velocities 
over  the  airfoil.  Since  the  theories  are  restricted  to  small  stream  shear 
(an  order  of  magnitude  smaller  than  the  experimental  conditions),  it  was 
not  possible  to  calculate  the  distributions  under  stream  conditions  duplica¬ 
ting  those  in  the  experiments.  However,  the  calculations  showed  two 
effects  tending  to  counterbalance  each  other.  The  zero-thickness  shear 
effect  is  to  cause  a  "saddle-backed"  distribution,  while  the  first-order 
shear-thickness  interaction  tends  to  cause  an  overall  decrease  in  velocity 
on  the  airfoil  under  conditions  similar  to  those  in  which  delayed  stall  was 
observed.  Since  the  former  is  larger  than  the  latter  effect,  within  the 
limits  of  the  theory,  stall  should  be  promoted. 

It  is  noted  that  parametric  dependence  of  the  destalling  effect  deduced 
from  the  experimental  results  (that  the  product  of  the  local  shear  and  the 
derivative  of  the  local  shear  be  negative)  is  of  higher  order  than  covered 
by  either  theory,  and  it  is  concluded  that  a  theory  valid  for  large  values  of 
stream  shear  is  required  to  check  the  experimental  results. 
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With  this  qualitative  understanding  of  the  cause  of  the  destalling  phen¬ 
omenon,  it  is  possible  to  speculate  on  a  number  of  practical  applications  that 
can  be  used  to  augment  wing  maximum  lift.  The  most  obvious  application  is 
to  exploit  the  shear  inherent  in  propeller-driven  aircraft  by  properly  orient¬ 
ing  the  wing  with  respect  to  the  slipstream  axis,  that  is,  the  wing  somewhat 
below  the  slipstream  axis.  This  application  might  be  particularly  effective 
in  deflected  slipstream  configurations  in  which  the  wing  is  almost  entirely 

immersed  in  a  shear  flow.  The  practicability  of  this  scheme  is  qualitatively 

1  2 

demonstrated  by  Brenckman's  experimental  results  and  by  more  recent 
CAL.  experimental  results  Both  show  that  the  destalling  phenomenon  is 

realized  with  axisymmetric  slipstreams. 

Another  potential  means  for  utilizing  the  destalling  phenomenon  asso¬ 
ciated  with  shear  flows  is  to  artificially  generate  a  shear  flow  ahead  of  a 
wing.  Previous  experimental  research  disclosed  that  stall  is  essentially  un¬ 
affected  in  two-dimensional  flow  by  the  proximity  of  the  slipstream  boundary. 
This  suggests  that  a  fairly  thin  shear  flow,  perhaps  comparable  in  thickness 
to  the  airfoil  thickness,  might  suffice  to  delay  stall.  Since  this  is  of  interest 
only  at  very  low  speeds,  the  shear  could  be  generated  with  nonuniform 
screens,  as  in  the  experiments.  There  would  be  a  drag  penalty  since  losses 
are  being  introduced,  but  the  losses  are  purely  viscous  and  vary  as  the 
square  of  the  velocity,  so  that  they  might  be  made  small.  Alternatively,  the 
shear  could  be  generated  with  air  jets  spanning  the  wing.  Again,  these 
would  be  employed  only  at  low  speeds  so  that  the  required  power  might  not 
be  prohibitive. 
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APPENDIX 


THE  THICKNESS-SHEAR  INTERACTION  IN  NONUNIFORM  SHEAR  FLOW 


A  potential  explanation  of  the  destalling  phenomenon  observed  in 

the  experiments  with  an  airfoil  in  nonuniform  shear  flow  is  that  the 

stream  shear  causes  alterations  in  the  airfoil  pressure  distribution.  This 

4 

possibility  was  investigated  by  using  the  available  theory  to  calculate  air¬ 
foil  pressure  distributions  under  conditions  similar  to  those  in  the  experi¬ 
ments.  That  is,  the  distributions  were  calculated  for  the  airfoil  located 
above  and  below  the  stream  plane  of  symmetry,  where  the  experiments 
showed  stall  should  be  promoted  and  delayed,  respectively.  The  calcula¬ 
tions  showed  the  two  pressure  distributions  were  nearly  identical  and 
suggested  stall  should  occur  at  the  same  lift  regardless  of  location  in  the 
stream. 

4 

It  was  noted  that  the  theory  was  approximate  in  that  the  stream 
shear  and  airfoil  thickness  were  assumed  small,  and  products  of  the  air¬ 
foil  thickness  and  the  stream  shear  parameters  were  neglected.  The 
experimental  data  show  there  is  a  sizable  thickness-shear  interaction, 
and  it  was  suggested  that  by  retaining  these  terms  in  the  theoretical 
pressure  distribution,  the  destalling  effect  might  be  explained.  This 
led  to  the  following  extension  of  Jones'  theory  to  predict  the  thickness- 
shear  interaction. 

In  a  flow  with  shear,  the  stream  function  is  governed  by  the 
equation 


v>  =  f(y^) 


(i-i) 


where  -f  is  an  unknown  function.  If  we  take  the  distribution  of  velocities 


in  the  remote  stream  to  be 


U  =  -  V  + 


±  x_. 

2  H* 


(1-2) 


the  corresponding  stream  function  is 
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Introducing  a  perturbation  stream  function  to  describe  the  flow  devia¬ 

tion  due  to  the  airfoil,  we  find,  to  the  first  order  in  <jL  ,  that 


f, 


(1-4) 


Equation  (1-4)  is  to  be  solved  subject  to  the  boundary  condition 


=  constant  on  airfoil  surface.  (1-5) 


The  boundary  value  problem  specified  by  (1-4)  and  (1-5)  is  now 
attacked  in  the  usual  wayf  As  a  first  step,  the  equations  are  transformed 
to  a  coordinate  system. aligned  with  the  airfoil  by  use  of  the  relations 


2  =  'ti  oc  —  ^  oc 

Y  =  6  +-  -1^  oi 


(1-6) 


The  geometric  parameters  e  and  oc  are  respectively  the  distance 
above  the  axis  of  symmetry  of  the  incident  flow  and  the  angle  of  attack. 
The  second  step  is  to  specify  the  airfoil  as  that  mapping  in  the  ^i-'y^plane 
from  the  circle  5=0  by  the  conformal  transformation 

•  {hi-  icj 

Z  =  ^  h  cj  =  Z^-h  +  -'Y  H  (1-7) 

®  e 


or 


+  H.  +(H-2ic„)€  -  <=• 


(1-8) 


which  is  correct  to  the  first  order  in  and  •y'«  .  The  geometrical 

parameters  introduced  in  Eq.  (1-8)  have  the  following  significance:  -j/, 
and  -y-,  are,  of  course,  the  coordinates  of  the  circle  of  radius  H  which 
transforms  into  the  airfoil,  but  in  addition,  specifies  the  camber  of 

the  wing.  For  a  symmetrical  airfoil  y-f/H  =  0  and  the  thickness  is  a 
unique  function  of  '^o/ ^  >  which  for  small  values  of  is  approxi¬ 

mately  linear. 
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Under  the  transformations  (1-6)  and  (1-8),  the  field  equation  (1-4) 


becomes 


-  t  .31 

:  /U>ay3^-hS  AU>^<P 


(1-9) 


v;  (i,f)  > 


while  after  considerable  manipulation,  the  boundary  condition  reduces  to 


K  f  +  ^(o)  .3^ A.  <p)  =  C  =  constant 


(I-IO) 


The  boundary  value  problem  specified  by  (1-9)  and  (I- 10)  is  difficult 

since  Eq.  (1-9)  is,  as  it  stands,  inseparable.  E.  E.  Jones,  in  Ref.  4, 

suggested  the  neglect  of  the  terms  in  ^  ,  which  reduces  the  problem 

"  Z  0  « 

to  the  partial  differential  equation  of  Mathieu  .  It  is  probably  possible 
to  regard  Jones'  solution  as  the  zero-oider  approximation  and  proceed 
by  perturbation  theory,  obtaining  the  equations 


-S 


(I-ll) 


for  the  subsequent  stages  in  the  iteration.  Equation  (I- II)  is  certainly 
easier  than  (1-9).  but  still  pjeseiits  some  puzzling  features.  The  form  of 
Eq.  (I-ll),  together  with  the  orthonormal  character  of  the  Mathieu 
functions,  suggests  some  adaptation  of  the  perturbation  procedure  of 
quantum  mechanics. 

While  this  boundary  value  problem  must  be  solved  at  least  approxi¬ 
mately  if  complete  information  on  the  effect  of  the  terms  is  to  be 

obtained,  some  information  can  be  obtained  without. 


Actually,  the  complete  solution  of  the  extended  boundary  value 
problem,  though  useful,  may  not  be  essential  for  our  immediate  purpose. 
Both  the  lift  and  the  separation  point  depend  on  the  tangential  velocity 
at  the  surface  of  the  airfoil.  This  tangential  velocity  is 


L_/^  ^  ^ 

Aif)  \  31  3l) 


(1-12) 


where  |  is  completely  independent  of  ^  .  Now,  under  the 

transformations  (1-6)  and  (1-8),  the  incident  stream  function  becomes 


yp  =  -V 

O 


(1-13) 


whe  re 

9)  OC 


(1-14) 


=  e  -h  Z  HyCt^dJv  ^  yC^9  ^M/rx^OC  2  // .-i-ixxA.  ^  -^Uyyty  ^  -C^^'oC 


Consequently,  the  thickness  effect,  ^  ,  is  seen  to  enter  both  terms  of 

Eq.  (1-12).  E.  E.  Jones,  in  Ref.  4,  neglected  the  last  term  in  (1-13)  as 
being  of  second  order  in  small  quantities,  .11  we  are  to  include 

this  term  we  must  add  an  increment 


=  - 


A(P) 


6*7^ -If-) 


1*0 


(1-15) 


to  the  tangential  velocity 


u--- 


A(9) 


4 

M?) 


(1-16) 


computed  by  Jones.  A  tedious  calculation  shows  that  the  velocity  increment 
can  be  expressed  as 

/  ^ 


C^U  = 


/\(f; 


6^  +  4?  ^  (p  ^  9 


^■-1 


(1-17) 


wh  ere 
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Now,  as  Jones  remarks  in  Ref.  4,  the  circulation  around  the  airfoil 
is  just  2  TT  times  the  constant  term  or,  in  our  case,  Z  rr  ~  ^ 
and  this  circulation  must  be  adjusted  by  varying  an  arbitrary  factor  in  ac^ 
to  make  the  velocity  finite  at  the  trailing  edge  of  the  contour.  In  this 
way,  the  circulation  and,  hence,  the  lift  is  seen  to  receive  a  contribution 
from  •  As  a  matter  of  fact,  a  repetition  of  the  argument  in 

Ref.  4  leads  to  the  conclusion 
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(1-19) 


2(*, - (5j(bj 


The  significance  of  Eq,  (1-19)  is  most  readily  determined  by 
examining  the  case  of  zero  incidence  (  oc.  «  O  ).  According  to  Ref.  4 


a,(o)  =ct,(0)^  a^iO)  =  ai(0)  =  b^iO)  =  b,(0)  =  0 


b,(o)  =  -v6w 


while,  according  to  Eq.  (1-18) 


S,  A  ^ 


/s,(o)  =  0 


Thus,  the  lift  increment  at  zero  angle  of  attack  is 


(1-20) 


(1-21) 


(1-22) 
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Figure  i  THEORETICAL  MODEL  OF  UNIFORM  SHEAR 


F^gure  2  THEORETICAL  MODEL  OF  NON-UNIFORM  SHEAR 
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gure  5  TWO-DIMENSIONAL  AIRFOIL  MODEL  AND  FORCE  BALANCE 


3.2 


Figure  6  WING  SECTION  LIFT  IN  A  TWO-DIMENSIONAL  PROPELLER  SLIPSTREAM 


DRAG  CRITERION  FOR  APPROACHING  STALL 
X  PITCHING  MOMENT  CRITERION  FOR  APPROACHING  STALL 


Figure  9  CORRELATION  OF  MAXIMUM  LIFT  WITH  STEAM  PROPERTIES 
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Figure  13  AIRFOIL  VELOCITY  DISTRIBUTION  IN  UNIFORM  AND  NON-UNIFORM  FLOW 
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